Density of states modulation in the pseudogap state of high-T c superconductors 
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The density of states modulation recently observed by the scanning tunneling microscopic experi- 
ment in the pseudogap state of high-T c superconductors is explained by the pairing assisted particle 
transitions under perturbation of a periodic pairing interaction. In such a transition process, a par- 
ticle with momentum k firstly picks up another particle of inverse spin to produce a pair leaving a 
hole. Under the perturbation of periodic pairing interaction of modulation wave vector Q, the pair 
absorbs a momentum Q and then breaks into two single-particles: one propagates with momentum 
k + Q, and another one fills in the hole. The transition is significant at low energies within the 
pseudogap since where the pairing excitations most favorably survive. We calculate the Fourier 
component of the modulated density of states using two different models. Both of the theoretical 
results are consistent with the experiment. 

PACS numbers: 74.25.Jb,74.72.-h,74.50.+r 



The recent scanning tunneling microscopic (STM) 
experiments on underdoped Bi2Sr2CaCu20g+5 samples 
have detected a dispersionless spatial modulation close 
to the periodicity of 4a x 4a in the tunneling conductance 
in the pseudogap states [1]. This modulation seems to be 
different from that in the superconducting state, which 
exhibits energy dispersion [2-4]. The STM observations 
show a distinct enhancement of intensity in the modu- 
lated patterns within the pseudogap energy scales. A 
similar modulation has been observed in the vortex cores 
[5]. This phenomenon provides a new clue for under- 
standing the physics of the pseudogap states. 

The mechanism of modulated electronic structures 
has been studied by various models of electron order- 
ing [6-13]. The possible origin of the modulation in the 
superconducting state is attributed to the quantum inter- 
ferences of quasiparticle scattering by some kind of per- 
turbation [11]. For explaining the dispersionless modula- 
tion in the pseudogap states, Chen et al. have proposed 
a model of density wave of d-wave Cooper pairs without 
global phase coherence [13]. The Fourier component of 
the tunneling density of state (DOS) p(Q, E) is predicted 
to be approximately even with respect to the energy E. 
This result captures the main feature of the energy de- 
pendence of the experimental observations. They noted 
that in case of charge density wave (CDW) or poten- 
tial scattering, p(Q,E) is approximately antisymmetric 
with E —> —E. Very recently, Anderson has suggested a 
model of Wigner crystal of d-wave hole pairs embedded in 
a background of d-wave resonance valence bond of singlet 
electron pairs and argued this structure of holes might be 
also relevant with the anomalous i-doping phenomenon 
[14]. 

In this paper, we will propose a description that the 
modulation phenomenon in the pseudogap state is a con- 
sequence of pairing assisted single-particle transitions un- 
der the perturbation of the modulated pairing interac- 



tion. Instead of supposing a crystal of pairs, we consider 
the pseudogap state as a system of single-electrons coex- 
isting with the uncondensed pairs. Such a description for 
the pseudogap state is consistent with the existed mod- 
els for studying the superconducting phase transition and 
the transport properties [15-20]. By model calculation of 
the Fourier component of the modulated DOS that essen- 
tially equivalent to the modulated tunneling conductance 
at low temperatures, we will show the present model can 
reasonably explain the STM experimental result. 
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FIG. 1. First order self-energy of the Green's function. The 
dashed line with a cross denotes the pairing interaction with 
modulation wave number Q. The shaded parts represent the 
ladder diagrams. 

The process of pairing assisted single-particle transi- 
tion is illustrated by the self-energy of the first order per- 
turbation as in Fig. 1, where the dashed line with a cross 
denotes the modulated pairing interaction of modulation 
wave number Q, and shaded rectangles are the ladder di- 
agrams representing in essence the pair propagation. The 
quantities k = (k, z n ) and q = (q, Z m ) are the general- 
ized momenta with z n and Z m the imaginary Matsubara 
frequencies respectively of fermions and bosons, q + Q is 
understood as (q+ Q,Z m ). In the process, a particle of 
momentum k firstly picks up another particle of inverse 
spin with momentum q— k to make up a singlet pair of to- 
tal momentum q. Under the perturbation of the periodic 
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pairing interaction, the pair transits to a state of momen- 
tum q + Q. When it breaks, it emits two single-particles: 
one fills in the hole, and another one propagates with 
momentum k + Q. As a result, the single-particle jumps 
from k to k + Q through the assistance of the pair tran- 
sition. Since the pair excitations are pronouncedly only 
within the pseudogap, such a single-particle transition 
and thereby the DOS modulation are significant within 
the same energy scales. In the following, we firstly il- 
lustrate that the DOS modulation is approximately even 
with E — ► —E, and then in the later part of this paper, 
we present our model calculations. 

At temperature T close to the transition point T c , the 
pair propagator has strong peak at q — 0, which im- 
plies the pair excitations happen mostly at low energies 
[16,19,20]. Especially this peak diverges at T = T c . 
To qualitatively analyze the DOS modulation, we here 
firstly consider the extreme case that the contribution 
to the first order self-energy E^^fcjQ) predominantly 
comes from small regions around q = and q + Q = 
(at which the right and left pair propagators respec- 
tively in Fig. 1 have strong peaks). We then have ap- 
proximately llW(k,Q) oc [G(-fc) + G{-k - Q)]/(k,Q) 
with /(k, Q) an interference factor [see Eq. (11)]. The 
first order Green's function behaves like G^'(k,Q) oc 
G(k + Q) [G(-k) +G{-k- Q)}G(k)f(k, Q). The Fourier 
component of the modulated DOS is defined by 



P(Q,E) = - 




where N is the total number of sites in a lattice on which 
the electron system is modeled. By the approximation 
above, we have 

p(Q, E) cx £ Im[G fl (k + Q, E)G A (-k, -E)G R (k, E) 

k 

+ (k<-k + Q)]/(k,Q). 

(2) 

The k-summation comes predominately from the regions 
close to the quasiparticle poles where the Green's func- 
tion behaves like 

G R (k,E) = [G A (k, E)]* « - ° k , (3) 

E - Ek + 27 k 

where E^ and 7k are respectively the energy and width of 
a quasiparticle, and ak is the residue at the pole. At low 
energy E <~ 0, only those states of E^ w -Ek+Q ~ E in 
some regions close to the Fermi surface are important to 
the k-summation. For example, Q = (tt/2,0) (in unit of 
lattice constant a = 1), these regions are near to (0, ±tt). 
When E is changed to —E, the corresponding states are 
changed to another side of the Fermi surface. Under the 
change G R (k,E) » -G A (k\-E) with E k , = -£ k , a 
minus sign comes from the three Green's functions in the 



square bracket in Eq. (2). This minus sign is just can- 
celed by the imaginary part of the complex conjugate. 
On the other hand, since /(k', Q) should not be differ- 
ent too much from /(k, Q) and the unperturbed densities 
of those important states of respectively k and k! are ap- 
proximately the same, we thus have p(Q, E) w p(Q, —E). 
From the discussion, we can draw an outline for p(Q, E) 
at small E: there is an approximately symmetric peak at 
E = with width of the pseudogap energy scale. 

To present our quantitative calculation, we start with 
the perturbation of the modulated pairing interaction, 

ij 

with Vij = Vo(|i"i — rjl) exp[Q • (n + r j ) / 2] and n ia the 
density operator of electron with spin-a at site i. By first 
order perturbation, since each Fourier component of the 
DOS modulation can be analyzed separately, we consider 
a single modulation wave vector and set Q = (n/2,0) in 
the present calculation. We will consider the singlet d- 
wave pairing of electrons at the nearest-neighbor sites. It 
is then convenient to write the perturbation Hamiltonian 
in k-space, 

H> = ^a? Y^^\+QPi + PqPq+Q)> ( 5 ) 
q 

where vq is the interaction strength of d-wave channel, 
and p q = X) 7 ?kCk+ q /2T c -k+q/2i is the d-wave pair op- 

k 

erator, with c ka (cu a ) the creation (annihilation) op- 
erator for electrons with momentum-k and spin-a and 
?7k = cos k x — cos k y . 

For the unperturbed electron system, we firstly con- 
sider a tight-binding model with d-wave channel interac- 
tion [16,17,19]. The Hamiltonian is given by 

ka q 

where £k = — 2t (cos k x + cos k y ) — 2t z cos k z — fj, with \i the 
chemical potential. For taking into account of the effect 
of short-range strong Coulomb repulsion, the hopping in- 
tegrals t and t z are considered as to be proportional to 
the doping concentration S, e.g., t = toS with to a con- 
stant. For the quasi-two-dimensional system, t z /t = 0.01 
is supposed. The weak interlayer coupling ensures a fi- 
nite transition temperature [16,19]. In the present calcu- 
lation, we take v/t = 0.2 and 6 — 0.125. Throughout 
this paper, we use the units in which h — ks = 1. 

The single-particle Green's function for the unper- 
turbed pseudogap state reads 

G(k,z n ) = [z„-a-S(k,z„)]- 1 (7) 

with E(k, z n ) the self-energy. For brevity, we will occa- 
sionally use the generalized momentum defined earlier, 
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e.g., G(k, z n ) is simply written as G(k). Within the 
model given by Eq. (6), the most important effect on 
the single-particle comes from the pair excitations. In 
the self-energy, we take into account this pairing effect 
through the ladder-diagram that is a conserving approx- 
imation. The self-energy is approximated as [19] 



£(fc) 



T 

N 



q/2 G(<? - k)P(q), 



(8) 



where P(q) is the summation result of the ladder- 
diagram, and is given by 



P{Q) 



v 2 n(q) 



(9) 



l + vU(q 

U ^ = ~Jf E vlG(k+,z n )G(k-,Z m - zn) (10) 

kn 

with k ± = k ± q/2. With the computation algorithms 
recently developed in Ref. [20], we have obtained a self- 
consistent solution to the Green's function. 

With the result for the unperturbed Green's function, 
the first order self-energy T,^\k, Q) can be obtained im- 
mediately. For the convenience of numerical calculation, 
we express T,^'(k, Q) in the following form by denoting 
in Fig. 1 the input and output momenta with k — Q/2 
and k + Q/2, respectively, 



(ii) 



with /(k,q,Q) = J7 k+Q/4 _ q/2 ?7 k _Q / 4_ q/ 2. By noting 
that J7k-q/2 can be factorized as ??k-q/2 = "0^ (k)t^(q) 
with V'^(k) = (cos k x , sin k x , — cos k y , — sin k y ) and 
(/?t(q) = [cos(g x /2),sin(g' a: /2),cos(g'j,/2),sin(gj,/2)], the 
q-summation in Eq. (11) can be manipulated by fast 
Fourier transform. The expression for the first order 
Green's function is given by 



Though the calculated p(Q, E) can reasonably reflect 
the feature of the experimental observation, the phe- 
nomenological model seems to be too simple. The model 
cannot correctly treat the short-range antiferromagnetic 
correlation between the electrons. Because of this reason, 
we have considered a more realistic model, the quasi-two 
dimensional Hubbard model. However, we will see that 
the results of the two models have the common feature 
as shown in Fig. 2; the qualitative behavior of the DOS 
modulation predicted by the process shown in Fig. 1 
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FIG. 2. p(Q,E) as function of E at various temperatures. 
($/&) depend ence on the perturbation strength vo is eliminated 
by dividing the constant 7 = vo/v. 

For the unperturbed state of the Hubbard model, we 
have recently obtained self-consistent solution using the 
charge, spin and pairing-fluctuation exchange approxi- 
mation [20]. By this approximation, besides the pair ex- 
citations are taken into account by the ladder-diagram, 
the effect of short-range Coulomb repulsion is treated by 
the standard charge and spin-fluctuation exchange. In 
case of the Hubbard model, the functions P(q ± Q/2) 
appeared in Eq. (11) arc modified to be P(q ± Q/2) = 



G (1) (fc,Q) = G(fc + Q/2)£ (1) (fc,Q)G(fc-Q/2). (12) v 2 U (q ± Q/2)/[l + vIL(q ± Q/2)} with 



By firstly integrating G^'(k, Q) with respect to the mo- 
mentum k and then using the Pade approximation for 
analytical continuation of the result from the imaginary 
Matsubara frequency to the real energy [21], we then ob- 
tain the result for p(Q, E). 

Shown in Fig. 2 are results for p(Q, E) at temper- 
atures T/t = 0.0175, 0.025, and 0.0375 all above the 
transition point. [The values of — vn(0) are 0.925, 0.827, 
and 0.700, respectively. The transition point corresponds 
to — vIl(Q) = 1.] As expected, there is a central peak in 
p(Q, E) nearly symmetric about E — 0. This is con- 
sistent with the experiment [1]. Out of the pseudogap 
region, p(Q, E) is negative, which should be consistent 
with a sum rule [10]. 



no(«) 



N 
T 



E%- q /2^(fc)G(fc)G(9-fc) (13) 



n (9) = -ivE^ 2 ( fc ) G ( fc ) G (9- fc ) 



(14) 



where 4>(k) is the pair wave function of ci-wave symme- 
try. Under a normalization condition for (j){k), the cou- 
pling strength v in this case is obtained as v — A/2 with 
A the eigenvalue of the Eliashberg equation [20]. Cor- 
respondingly, the interference factor /(k, q, Q) in Eq. 
(11) should be changed to /(k,q, Q) = <f){k + Q/4 - 
q/2)4>{k — Q/4 — q/2). However, with such an inter- 
ference factor, it is nearly impossible to get a numer- 
ical result for S^)(fc, Q); the amount of the computa- 
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tion is tremendous. We then make an approximation by 
factorizing <j>(k ± Q/4 - q/2) oc A(z„)77 k± Q/ 4 _ q / 2 with 
A(z„) = 0(kM, z n )/2 and kjvi = (tt,0). That means the 
wave function of a pair has approximately the nearest- 
neighbor cZ-wave pairing structure and its frequency de- 
pendence is separated by the factor A(z„). The pro- 
portional coefficient is determined by the normalization 
condition 

^^A 2 (z„)^ G(fc)G( _ fc) = A (15) 

k 
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FIG. 3. p(Q, E) as function of E for the Hubbard model. 
The constant 7 is given by vo/v p . 

In Fig. 3, we depict the result for p(Q, £7) for the 
Hubbard model at U = 5t, 5 = 0.125 and T = 0.0123*, 
where A = 0.982 and v p — 15.5 are obtained. The main 
feature of this result is approximately the same as in Fig. 
2; both of them show a central peak at E = 0. From 
Fig. 3, the width of the peak is seen about 0.2* that is 
overall the same as the pseudogap width found in the 
unperturbed DOS [20]. This is again consistent with the 
experiment. The central peak seems less symmetric than 
expected. This is because the unperturbed DOS of the 
Hubbard model by the approximation is not symmetric. 
Near the Fermi energy, the unperturbed DOS at E < is 
larger than that at E > 0, which means the summation 
in Eq. (1) over the states below the Fermi surface is 
larger and results in a shift of the maximum of p(Q, E) 
toward to a slightly negative energy. 

In summary we have shown that the DOS modulation 
observed in the pseudogap state of high-temperature su- 
perconductors can be explained by the pairing assisted 
particle transitions under the perturbation of modulated 
pairing interaction. Such a transition process is illus- 
trated in Fig. 1. The calculated Fourier component 
p(Q, E) as a function of energy has a central peak at 
E = with a width of pseudogap energy scale in consis- 



tent with the experiment. 
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